Structure factor of low-energy spin excitations in a S = 1/2 kagome antiferromagnet 
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The ground state of the 5 = 1/2 Heisenberg antiferomagnet on kagome can be viewed as a col- 
lection of fermionic spinons bound into small, heavy singlet pairs. Low-energy magnetic excitations 
in this system correspond to breaking the pairs into individual spinons. We calculate the structure 
factor for inelastic neutron scattering from independent spinon pairs. 



I. INTRODUCTION 

The spin- 1/2 Heisenberg antiferromagnet on kagome 
(Fig. Q]) has attracted attention of theorists for more 
than two decades^* 2 -. Synthesis of new magnetic ma- 
terials with Cu 2+ spins forming precisely such a lat- 
tice has generated renewed interest in this systemj^ 
Strong frustration, induced by the non-bipartite geom- 
etry, and strong quantum fluctuations suppress long- 
range magnetic order and open intriguing possibilities 
of an exotic quantum-disordered ground state and un- 
usual magnetic excitations. Numerical calculations 5,6 
suggest that the ground state of the system is a total 
spin singlet and that magnetic excitations are gapped. 
Among the proposed ground states of this system are a 
valence-bond crystal (VBC}£~— breaking lattice symme- 
tries or a spin liquid of some sorti^Sr— Recent synthe- 
sis and experimental characterization of herbertsmithite 
ZnCu 3 (OH) 3 Cl 2 r^ where Cu 2+ ions carry S = 1/2 
and form a perfect kagome lattice, provides additional 
motivation for theoretical studies of this model. 

The simplest model takes into account exchange in- 
teractions between nearest neighbors represented by the 




FIG. 1: A dimer covering of kagome. Defect triangles (D) 
lack dimers. 



Hamiltonian 

H = jJ2 S i- S j> (!) 

<ij> 

Recasting the Hamiltonian as JJ2a ^a/^, where Sa is 
the net spin of triangle A, suggests a plausible route to 
constructing the ground state of the system: minimize 
the total spin of every triangle, Sa = 1/2, by locking two 
of its three spins in a S = bond. Unfortunately, this 
program fails, as a simple counting argument showsi A 
kagome with N triangles contains 37V/2 spins and thus no 
more than 3N/4 singlet bonds can be formed, leaving at 
least one in four triangles without a singlet bond (Fig. [I}. 
As a result, a short-range valence-bond state is not a 
ground state of the Hamiltonian (QJ. 

Although the naive approach to constructing a valence- 
bond ground state fails, it nonetheless provides a useful 
starting point by forcing us to regard most of the lattice 
(triangles with a singlet bond) as the vacuum contain- 
ing relatively dilute dynamical objects (defect triangles). 
Elucidating the properties of an isolated defect triangle is 
the next logical step towards understanding the physics 
of the model. Work in this direction was carried out early 
on by Elser and Zengj 2 -^ who considered an isolated de- 
fect triangle on the Husimi cactu a 25 i 26 , a Bethe lattice 
of corner-sharing triangles. They found that quantum 
fluctuations are strongly localized in the vicinity of a de- 
fect. By treating the kagome antiferromagnet as a dilute 
ensemble of fluctuating defects, Elser and Zeng obtained 
an estimate of its ground-state energy in excellent agree- 
ment with numerical diagonalization. 

More recently, we pointed out^ 7 . that defect triangles 
are composite objects: they are bound states of two 
quasiparticles with 5* = 1/2. These quasiparticles bear 
strong resemblance to spinons of the A, or sawtooth, 
chai n 28 i 29 . The spinons come in two flavors: kinks and 
antikinks. Kinks are localized and have zero excitation 
energy, while antikinks are mobile with the lowest en- 
ergy of 0.218J. A defect triangle on the Husimi cac- 
tus is a tightly bound state of two antikinks with to- 
tal S = 0. Lowest-energy spin excitations correspond 
to breaking up the pair into two mobile antikinks. The 
spin gap is set by the binding energy of the pair, 0.06J. 
This value, obtained for a single defect on a cactus, 
agrees with the spin-gap estimates for kagome based 
on a series expansion 7 (0.08J ± 0.02J) and DMRG& 
(0.055J± 0.005 J). This agreement suggests that view- 
ing the kagome antiferromagnet as a dilute ensemble of 
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spinon pairs is a useful point of departure for further ex- 
plorations of this model. 

We would like to stress that antikink spinons are not el- 
ementary excitations of the kagome antiferromagnet but 
rather its building blocks. (In the same way, quarks are 
building blocks of a baryon but not its elementary exci- 
tations.) One in four triangles carrying a spinon pair 
translates into 1/3 of an antikink per site. We have 
found that antikinks exhibit fermionic statistics. The 
minus sign upon an exchange of two antikinks comes 
from the "Berry phase" associated with the adiabatic 
motion of singlet bonds. While this picture of fermionic 
spinons resembles some of the previous proposals inspired 
by fermionic large- N expansions (e.g., Hastings^), there 
are important differences. We find that fermions expe- 
rience a strong exchange-mediated attraction in the sin- 
glet channel, which causes them to form tightly bound 
pairs, thus invalidating the picture of a fermion sea with 
a Fermi surface or Dirac points. We also find that the 
emergent compact U(l) gauge field manifests itself not as 
a background magnetic flux, but as a quantized electric 
field of unit strength whose presence strongly constrains 
the motion of antikinks carrying a U(l) charge of +2 and 
even more strongly affects antikink pairs (charge +4) . It 
remains to be seen whether further progress can be made 
in the problem of spinons interacting with one another 
and the gauge field, a many-body problem with strong 
interactions. 

In the present work, we calculate the dynamical struc- 
ture factor of low-energy spin excitations in a kagome an- 
tiferromagnet observable by inelastic neutron scattering. 
The physical process responsible for the lowest-energy 
magnetic scattering is the breaking up of a S — an- 
tikink pair into two antikinks with parallel spins that 
subsequently move away from each other. If the ground 
state of the system is a valence-bond crystal the mov- 
ing spinons disturb the preferred valence-bond arrange- 
ment. The resulting energy increase leads to spinon 
confinement. We neglect this effect because numerical 
diagonalization— and series expansion^ indicate that the 
energy differences between various valence-bond config- 
urations are very small, on the order of 10 _3 J per site. 
We therefore expect that the confinement length is long 
and that its effects can be neglected in the first approxi- 
mation. 

We also neglect the influence of antikink pairs on one 
another. Although the liberated antikinks may run into 
other pairs present on the lattice, this is not a severe 
problem. The structure factor is determined by the over- 
lap of the initial and final wavefunctions of the antikink 
pair. Because the initial state is well localized, the result 
of the calculation is not sensitive to the long-distance 
behavior of the final state. We compute the dynamical 
spin correlations in real space in the presence of a single 
antikink pair on the Husimi cactus. The resulting cor- 
relation function decays quickly as we move away from 
the defect triangle. Assuming that the spin correlations 
have similarly local nature on kagome proper, we trans- 




FIG. 2: A correspondence between spinon trails on the Husimi 
cactus (a) and kagome (b). Spins along the three shaded trails 
are labeled San. Here a = x,y,z denotes a trail and n — 
1, 2, 3 . . . enumerates the spins along it. The correspondence 
breaks down for n > 5 because trails on kagome may begin to 
overlap. Further complications are brought by the presence 
of other defect triangles (labeled D) next to the trails. 



late the obtained spin correlations to kagome using the 
correspondence depicted in Fig. Finally, after a spa- 
tial Fourier transform, we obtain the structure factor in 
k space at the edge of the spin gap. 

The paper is organized as follows. In section [TlJ we 
review the calculations of the wave functions of two an- 
tikinks with total spin 5 = and 1 on the Husimi cac- 
tus. Section lnll describes the calculation of the dynamical 
structure factor. We discuss the results in section ITVl 
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FIG. 3: Valence-bond states of the Husimi cactus. The state 
with a single defect triangle (a) evolves into states with two 
antikinks (b-d) traveling through the cactus. The dashed line 
indicates that the two antikinks have total spin 0. 



II. CALCULATIONS OF WAVE FUNCTIONS 

In order to study a defect triangle in isolation, Elser 
and ZengS examined kagome lattices in two-dimensional 
spaces of constant curvature where hexagonal loops of 
the familiar kagome are replaced with loops of length 
L. Lattices with L < 6 are finite and can be embed- 
ded in a two-dimensional sphere^ whereas lattices with 
L > 6 are infinite and live in two-dimensional hyperbolic 
spaces^ The L = oo structure is a tree of corner-sharing 
triangles known as the Husimi cactu s 25 ' 26 (Fig. [2]). It 
has just the right ratio between the numbers of sites and 
links to permit the construction of a static valence-bond 
ground state and to study a defect triangle in isolation. 

A single defect triangle turns out to be a bound state 
of two antikink spinons with total spin 0. The bound 
state is formed because spinons with antiparallel spins 
experience exchange-mediated attraction. In contrast, 
two antikinks with total spin 1 repel and thus do not 
form a bound state. 



A. Two antikinks with 5 = 

1. Hilbert space 

We begin by characterizing the Hilbert space to which 
the defect triangle state belongs. Upon applying the ex- 
change Hamiltonian on the defect triangle state once, the 
defect triangle breaks into two defects of a new type that 
connected by a long range singlet: for example, if we ap- 
ply left bond of the center trianlge in part a of Fig, the 
resulting state is shown in part b of Fig. These defects are 



close analogs of antikinks found in the one-dimensional 
sawtooth, or A chain . 28 ' 29 They are domain walls sepa- 
rating the two distinct ground states of the A chain. The 
antikinks carry spin 1/2 and are thus spinons. Under the 
action of the exchange Hamiltonian, the two antikinks 
can move along three one-dimensional paths that meet at 
location of the defect triangle. We will identify the three 
trails as x y and z (Fig. 2]) . The two antikinks can never 
be on the same trail. A generic state can therefore be 
written as \x,y, z), where x, y, and z are integers whose 
product vanishes. For example, in the state |2,3,0), one 
of the antikinks is on the second triangle of the x trail, 
while the other is on the third triangle of the y trail. 
Note that states |0,0, 1) |0, 1, 0) or 1 1,0,0) are identical 
to 1 0,0,0) (the original defect triangle), so they can be 
omitted altogether. It is the mathematical equivalence 
of the fact that any of the three singlets adjacent to the 
defect triangle can be viewed as a bound state of two 
anti-kinks. 

The resulting Hilbert space is denoted A°- Generally 
is the Hilbert space of n antikinks with total spin S. 
The Hilbert space A\ is illustrated in Fig. SJ 

Two basis states |r) = \x, y, z) and |r') = \x' , y' , z') are 
not orthogonal to one another because they are not eigen- 
states of the same Hermitian operator. Their overlap can 
be computed by following the standard recipe^ 



r r = 



(2) 



where the product is over loops formed by superimpos- 
ing if the dimer coverings of states |r) and |r'), L c is the 
length of loop c, and e c = ±1 is a Zi phase factor depen- 
dent on the sign convention for spin singlets. The S = 
state of two spins on sites i and j is antisymmetric under 
exchange: \(i,j)) = (\ Uj) ~ I UW)/^= HO'.*))- To 
remove the ambiguity, the state \{i,j)) is shown as an 
arrow pointing from i to j. The singlet phases need to 
be specified on the central triangle and along the three 
trails; all the other singlets are localized and their 
phases are irrelevant. Because the three trails can be re- 
garded as sawtooth chains, we choose the following phase 
convention.— A singlet points from the base of a triangle 
on the sawtooth chain to its vertex. In states \x,y,0), 
|0, y, z) and \x, 0, z), the long-range singlet has the same 
direction with the singlet on the center triangle, as shown 
in Fig. |H In states |x,0, 0), |0,y, 0) and |0, 0, z), the long- 
range singlet points from the central triangle to the trian- 
gle on one of the trails. Overlaps of various valence-bond 
states have the following values: 



(x,y,0\x',y',0) = 2 -l*-*'Hir-«r'l j 

(x,o,o\o,v,z) = 2 2 ~*-y-*, 

(x,y,0\x',0,z') = 2-1— '!-»-«', 

(x,t/,0|0,0,0) = 2 1 ~*-\ 

(x,0,0|0,0,0) = 2 1 -*. 

The rest of the overlaps can be obtained by using per 
mutations of x, y, and z. 



(3a) 
(3b) 
(3c) 
(3d) 
(3e) 
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FIG. 4: Hilbert space A° of two antikinks with S = 0. (a) Two 
antikinks can move among three trails of the cactus labeled 
as x, y, and z trails, (b) The lattice of \x,y,z) states with 
xyz — resembles the corner of a cube. 



With two antikinks in the three trails, at least one of 
the "coordinates" x, y, z is zero. For that reason, the 
motion of two antikinks can be mapped onto the motion 
of a single particle on three faces of an infinite cube, Fig. 
2] (b). This mapping reveals a simple rule that can be 
used to determine the overlap in all cases described in 
Eq. ([3]). Let n(r,r') be the length of the shortest path 
(expressed as the number of links) that connects points r 
and r' in Fig. @] (b) . The overlap between the two states 
is 



2. Orthonorrnal basts 

The basis states \x, y, z) can be orthogonalizcd by per- 
forming a simple rotation. We will denote the orthogo- 
nalized basis states as \x,y,z)° and the orthonormalized 
states as \x, y, z) on . In the orthogonalization procedure, 
the defect triangle state remains the same: 

10,0,0)° = 10,0,0). (5a) 
Its nearest neighbors are transformed as follows: 



|l,l,0)° = |l,l,0)-i|0,0,0). 



(5b) 



The resulting states |0, 1,1)°, |1,0,1)°, and |1,1,0)° are 
orthogonal to |0, 0, 0)° and to one other. 

The transformation for states along cube edges is: 

|M,0)° = |a;,0,0)-i|a: -1,0,0). (5c) 

The resulting state \x, 0,0)° is orthogonal to any state 
|0, y, z) with no antikink in the x trail. It is also orthog- 
onal to any state \x' , y, 0) with an antikink closer to the 
origin, x 1 < x. From that it follows that edge states 
\x, 0,0)°, \0,y, 0)°, and |0,0, z)° are orthogonal to one 
another and to the four states with antikinks near the 
origin. 

Finally, states \x, y, 0) with antikinks away from the 
origin (x,y > 1) are transformed by combining transla- 
tions in two directions: 



\x,y,0}° = \x,y,0)-±\x-l,y,0) 



(5d) 



~\x,v-l,0) + ^\x-l,v- 1,0). 



The orthonorrnal basis {|x, y, 0) on } is now easily ob- 
tained by normalization: 



0,0,0) on = 


|0,0,0)°, 


(6a) 


1,1,0)™ = 


^im,o)°, 


(6b) 


x,0,0) on = 


4m, o)°, 


(6c) 


x,y,0r = 




(6d) 



(r|r') = 2~"^ r '). 



(4) 



3. Effective Hamiltonian and the spectrum 
In the bulk, the Hamiltonian matrix is very simple: 
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H\x,y,Oy 



l,y,0) t 



J C.J 

-\ x ,y-l,0)°»+—\x,y,0)°° (7a) 



The Hamiltonian can be loosely understood as follows: two particles(antikinks) hop on the three one dimensional 
trails with — J/2 as their hopping amplitude and an attractive short-range potential. (We measure the energy from 
its value in the ground state.) 

The action of the Hamiltonian on the state with a defect triangle is as follows: 



i?|0,0,Or--^|l,l,0) on -^|l,0,l) on -^|0,l,l) on + ^|0,0,0) c 



(7b) 



(Note the reduced diagonal term.) 

For a state with one of the antikinks next to the defect triangle and the other farther away, \x, 1, 0) on with x > 1, 
y = 1, we have 



h\x, i, or = ~\x + 1, i, or - ^\x - 1, i, or - 2, or - ^-\ x . 0. ■ 



5J, 



1,0) c 



J 



4 2 4 

For a state with one antikink on the original defect triangle and the other somewhere in trail x, \x, 0, 0) c 



la:, 0,1) 



#|a:, 0, 0) on =a- + 1,0,0) 



3 -\x 1.0.0' 



x,l,0) c 



x,0, l) on + 2J|x,0,0) c 



(7c) 



(7d) 



Other special cases are considered in the supplemental 
material of Ref. H3- 

To determine the energy spectrum, we truncated the 
cactus at a finite radius R and diagonalized the Hamil- 
tonian numerically. The spectrum for a singlet pair of 
antikinks consists of a two-particle continuum starting at 
E = J/2 and a non-degenerate bound state at E — 0.44 J, 
i.e. 0.06 J below the bottom of the continuum. The di- 
ameter of the bound state, obtained by fitting its energy 
E(R) to the form E + SE exp(-2R/£), is £ = 2.8 lattice 
spacings. The probability to find both antikinks on the 
same triangle is 0.72, which implies a tight bound state. 



the overlap matrix elements between the states in A\: 



Two antikinks with S 



The energy spectrum of two antikinks in a S — 1 state 
was computed in the same fashion. We mapped out the 
Hilbert space A\ for a triplet pair by breaking one of 
the three singlets adjacent to the defect triangle at the 
center and letting the resulting spinons propagate along 
the three trails. The basis states, also be denoted \x,y,z), 
are again non-orthogonal. This time, there is no |0, 0,0) 
state, while |1,0, 0), |0, 1,0) and |0, 0,1) states do exist 
and are physically distinct. Here are some examples of 



(xi,yi,Q|0, 2/2,0) 

(xi, yi, 0|a;2,0,0) 

(aii, 2/1, 0|x 2 , 2/2)0) 
(xi,2/i,0|0,0, z 2 ) 

(xi, 2/i, 0|0, j/2, z 2 ) 



xi + \yi-V2\ 



yi+\xi-x 2 \ 



\xi-x 2 \ + \yi-y2\ 



Xl+\yi-V2\+Z2 



(8a) 

(8b) 

(8c) 
(8d) 
(8e) 



Other cases can be obtained by making use of the per- 
mutational symmetry. 

A similar procedure can be followed to orthogonalize 
most of the states. However, we did not find an obvious 
way to orthogonalize the states with spinons next to the 
original defect triangle, namely 1 1,0,0), 1 0,1,0), |0, 0, 1), 
1 1, 1, 0) , 1 1 , 0, 1) and |0, 1,1). We orthogonalize the rest 
of the basis states as exemplified by 



|x,l,0)°Hx,l,0) + i|x,0,0) 



(9) 



~2 lX ' 



1,1,0) --Ix- 1,0,0). 



Because some of the basis states are not orthogonal, the 
energy spectrum is obtained by solving the generalized 
eigenvalue problem: Hip = EOip where H is the Hamil- 
tonian matrix, O is the overlap matrix and tp is an eigen- 
vector. 
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FIG. 5: Hilbert space A\ of two antikinks with S = 1. (a) 
Minus signs denote bonds with the inverted sign of the hop- 
ping amplitude, t > 0. (b) Making a loop around the center 
changes sign of the wave function. 



The low-energy effective Hamiltonian is similar to that 
of the S = sector, with two important differences. 
First, the hopping amplitude between two states with 
a negative overlap is positive instead of negative. Most 
of these "wrong" signs can be corrected by a gauge trans- 
formation. However, the total phase accumulated on a 
loop containing the center in Fig. [5] adds up to a flux 



it, which reflects the fermionic nature of the antikinks. 
Second, in the singlet case the "potential" energy of the 
two antikinks is lowered by 3 J/4 when they both reside 
on the central triangle. In contrast, antikinks with par- 
allel spins repel each other: the energy of the system 
is raised by J/4 when they share the central triangle. 
Consequently, the spectrum of a triplet pair of antikinks 
consists of a two particle continuum starting at E = J/2 
with no bound states. The ground state in the S = 1 
sector is doubly degenerate. Both of these ground states 
have a line of nodes. The symmetry and degeneracy of 
the ground states are consistent with the Fermi statistics 
of antikinks. 



III. DYNAMICAL STRUCTURE FACTOR 
CALCULATION 

The real-space dynamical structure factor is defined as 

S(u, R,R') = J2 S(E f - Ei - hu)(i\Sk\f)(f\S+,\i). 
f 

(10) 

Inelastic neutron scattering directly measures its Fourier 
transform, 



5jv(w,q) 



/ 



5(Ef -Ei- Hlo) 



E 

R 



In both the initial state |i) and in final states |/), the mo- 
tion of the antikinks is restricted to three one-dimensional 
paths. Singlet bonds along these trails shift as the 
spinons move past them. In contrast, singlet bonds off 
the trails remain stationary. If R is located off a spinon 
trail, states S^ji) and |/) are orthogonal because the for- 
mer has a triplet on a bond involving site R, whereas the 
latter has a singlet there. For that reason, the struc- 
ture factor (TlU)) is nonzero only when both R and R' 
are on the spinon trails. Spins along the trail will be la- 
beled S an , where a = x, y, or z is the trail index and 
n = 1, 2, 3 ... is the position along the trail (Fig. [2]). R Q „ 
is the physical location of that spin on kagome. 

We evaluate the real-space structure factor (TTU|) on the 
cactus for the lowest energy transfer equal to the binding 
energy of an antikink pair. These dynamic correlations 
turn out to be sufficiently short-ranged to justify the ex- 
trapolation of the result to a pair of antikinks on kagome 
via the correspondence illustrated in Fig. [5J We assume 
that antikink pairs are randomly distributed over the lat- 
tice and therefore the net scattering intensity is the sum 
of scattering intensities of individual pairs (no interfer- 
ence). Because different pairs live in different valence- 
bond backgrounds, we average the structure factor over 
all possible configurations of the spinon trails. 

To simplify the calculations, we only approximate the 
singlet ground state as a superposition of four states 
|0,0, 0;s) (the defect triangle), |l,l,0;s), |0, 1,1; s) and 
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FIG. 6: The computed dynamical structure factor iV = 2, 3, 4, 
and 5 as well as for a few isolated dimers randomly distributed 
on the lattice. Wavenumbers are measured in units of 2n/a, 
where a is the lattice constant. 



1 1, 0, 1; s) (spinons hop to adjacent triangles). The trun- 
cated ground state has the overlap of 0.92 with the actual 
one. No truncation is done for the final triplet states, in 
which the spinons are delocalized. 
The form factor now reduces to 



f 



(12) 



N 



In this expression, we include spins S an located along 
the trails with n < N. In theory, one would like to 
take the limit N — > oo. In practice, we can only go to 
N = 5; for larger values of N trails may start overlapping 
on kagome. Fortunately, it can be seen in Fig. [6] that the 
structure factor changes little between N = 4 and N = 5, 
so we use the N = 5 result as our final answer. 

The matrix elements (f\S^ n \i) are calculated as fol- 
lows. The initial singlet ground state is a linear combi- 
nation of four states \x, y, z; s) from the singlet space A%. 
The final states are linear combinations of \x,y,z;t) in 



the triplet space A\. (As in previous sections, we use 
the shorthand r for the three "coordinates" of the an- 
tikink pair x,y and z. The physical coordinates of a site 
on kagome are denoted R Qn .) We first expand the initial 
and final states in the non-orthogonal singlet basis {|r; s)} 
and partially orthogonalized triplet basis {|r;t) on }: 



\i) = E c r( s )l r ; s )' 

r 

l/> = E^/WM)™ 



(13a) 
(13b) 



where / labels the two final states. After a computation 
of matrix elements 

M Q „(r,r')^ on (r;t|S+ n |r';s>, (14) 
listed in the Appendix, we obtain 

(f\S+ n \i) =E</( t ) M ««( r ' r ')cr(s). (15) 

r,r' 

and use it in Eq. (fT2")l . 

The dynamical structure factor is displayed in Fig. [5] 

IV. DISCUSSION 

We have computed the q dependence of the dynami- 
cal structure factor 5(q, uj) at the edge of the spin gap, 
A = 0.06J. Excitations responsible for this spectral 
weight are pairs of antikinks with parallel spins freed from 
their 5 = bound state. The wavevector dependence of 
the dynamical structure factor resembles that of isolated 
pairs of spins interacting via Heisenberg exchange: 



S{ui, q) = E 2I 1 ~ cos (q ' a)]5(ftw - J), 



(16) 



where a are separations of spins within a pair. For three 
dimer orientations, as on kagome, S(lj, q) is shown in 
Fig. [5] next to our result for the structure factor resulting 
from the breaking of antikink pairs. The resemblance is 
not surprising because the ground state of the Heisenberg 
model on kagome can be pictured as a collection of slowly 
resonating valence bonds. The structure factor of the in- 
elastic neutron scattering in powder herbertsmithite in- 
deed bears resemblance to that of isolated dimersj22 A 
similar q dependence was found for the instantaneous 
(frequency-integrated) structure factor by Laeuchli and 
Lhuillier— by using exact diagonalization. More recently, 
Singh 33 computed the dynamical structure factor due to 
another kind of excitations, viz. creation of kink- antikink 
pairs with a larger onset energy, A' = 0.25J, and found 
a similar intensity distribution. 

Our calculation is based on several assumptions. (1) 
We neglect the effects of virtual excitations in the form 
of kink-antikink pairs near a moving antikink. Our pre- 
vious study^i of a related one-dimensional system (the A 
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chai n 28 ' 29 ) showed that these excitations create a small 
renormalization of the antikink parameters but are oth- 
erwise harmless. (2) We assumed that the bound pairs 
of two antikinks with S = are sufficiently small so that 
their properties on a periodic kagome lattice and its tree- 
like analog are essentially the same. The small diameter 
of a pair, £ = 2.8 lattice spacings, provides assurance 
that this approximation is not unreasonable. (3) We ne- 
glected the dynamics of the pairs in the ground state of 
the Heisenberg model on kagome. This assumption is 
justified if the pairs move slowly on the time scale of the 
inverse spin gap. Some indications that this is so come 
from numerical work indicating a large density of sin- 
glet states at low energy^ Series expansion also indicates 
that the energy splittings between low-lying singlet states 
are very small, with energy differences as small as 10 -3 J 
per site. 34 (4) We assumed that pair positions exhibit 
no long-range order and thus their scattering amplitudes 
add incoherently. This may not be the case if the system 
indeed has valence-bond order, which requires a periodic 
arrangement of spinon pairs in the ground state. How- 
ever, the closeness of energy levels in the singlet sector 
suggests that the valence-bond crystal is fragile and can 
easily turn into a disordered solid under the influence 
of bond disorder or nonmagnetic impuritiesi 33 ' 35 (5) We 
ignored the possibility of interactions between adjacent 
pairs of antikinks. The dynamics of antikinks may be 
altered if their trails pass close to another defect triangle 
(Fig. [5]). This is a many-body problem that we hope to 
address in future work. 
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Appendix A: Matrix elements 

Before we give a list of matrix elements M an (r; r'), we 
need to specify a convention of labeling all the spins. The 
three one-dimensional trails on which the antikinks are 
moving can be viewed as three sawtooth chains originat- 
ing at the center triangle. We name a spin using the trail 
it's on and the distance between it and the center trian- 
gle. The three vertices of the center triangle are labeled 
S xX , S B i and S zi , Fig. [2] 

Since we write the initial state a linear combination 
of |0,0,0;s), |l,l,0;s), |0,l,l;s) and |l,0,l;s) and our 
convention for basis states respects the symmetry of per- 
mutations among "coordinates" x, y, and z, we only list 



nonzero matrix elements involving states |0, 0,0; s) and 
|l,l,0;s). For |0,0,0;s), 



(l,0,0;t|S+|0,0,0;s 
(0,l,0;t|5+|0,0,0;s 
(0,0,l;t|5+|0,0,0;s 

° n <n,0,0;t| 1 S+ 2B _ 1 |0 1 J 0;s 

oa <0,n,0;t| J S+ 2B _ 1 |0,0,0;s 

on (0,0,n;t|,S'+ 2 „_ 1 |0 ) 0,0; S 

with n > 2. Similarly for |1, 1,( 

on (n,l,0;t|5+ 2n _ 1 |l,l,0;s 
on <n,0,0;t|S+ 2B _ 1 |l,l,0;s 
oa <l,n,0;t|5+ 2n „ 1 |l,l,0;8 
™<0,n,0;t|S+ 2n „ 1 |l,l J 0;s 
° n (0,0 J n;t|S+ an _ 1 |l,l,0;s 



(1,0, 0;t 
(0,1, 0;t 
(1,0, l;t 
(0,1, l;t 
(1,0, 0;t 
(0,1, 0;t 
(l,0,l;t 
(0,1, l;t 
(1,1, 0;t 



|5+|l,l,0;s 
|5+|l,l,0;s 
|<S+i|l,l,0;s 
1^11,1, 0;s 
|5+|l,l,0;s 
|5+|l,l,0;s 
|5+|l,l,0;s 
|5+|l,l,0;s 
|5+|l,l,0;s 



1 

~ "71' 
i 

" "71' 

1 

= "71' 
= --(- 

= --(- 

---(- 
V6 \2 



n-1 



,(A1) 



V2G; 



71—1 



= -±a 



V6 V2 

7lG; ' 

-(- 

V6 V2 

-(- 

V6 V2 

1 



71—1 



2V2' 
1 



2>/2' 
1 



4V2' 
1 

1 



1 



2V2' 
1 



4V2' 
1 



2y/2' 
1 

"71' 



(A2) 
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